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ABSTRACT 

Recent proper motion studies of stars at the very center of the Galaxy 
strongly suggest that Sagittarius (Sgr) A*, the compact nonthermal radio 
source at the Galactic Center, is a 2.5x10^ Mq black hole. By means of 
near-simultaneous multi-wavelength Very Long Baseline Array measurements, 
we determine for the first time the intrinsic size and shape of Sgr A* to be 
72 Rsc(*) by < 20 Rsc(*), with the major axis oriented essentially north-south, 
where Rsc(*) (= 7.5 x 10^^ cm) is the Schwarzschild radius for a 2.5x10^ Mq 
black hole. Contrary to previous expectation that the intrinsic structure of 
Sgr A* is observable only at A < 1 mm, we can discern the intrinsic source size 
at A 7 mm because (1) the scattering size along the minor axis is half that along 
the major axis, and (2) the near simultaneous multi-wavelength mapping of 
Sgr A* with the same interferometer makes it possible to extrapolate precisely 
the minor axis scattering angle at A 7 mm. The intrinsic size and shape place 
direct constraints on the various emission models for Sgr A*. In particular, the 
advcction dominated accretion flow model may have to incorporate a radio jet 
in order to account for the structure of Sgr A*. 

Subject headings: galaxies: active — Galaxy: center — scattering 
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1. Introduction 

Sagittarius (Sgr) A*, the extremely compact non-thermal radio source at the Galactic 
Center, has been for many years considered the signpost of a massive black hole (e.g. 
Lynden-Bell and Rees 1971; Lo et al. 1985). Recent proper motion and radial velocity 
measurements of the stars in the immediate neighborhood of Sgr A* have provided very 
compeUing dynamical evidence for the existence of a compact dark mass of 2.5x10^ 
located within 0.015 pc (4.5x10^^ cm) of Sgr A*, supporting the hypothesis that Sgr A* is 
powered by a single massive black hole (Eckart and Genzel 1996; Eckart and Genzel 1997; 
Ghez et al. 1998). Since Sgr A* would be the nearest example by far of such a system, 
determining its intrinsic source structure would be very important for probing the region 
immediately surrounding the massive black hole. 

Up till now. Very Long Baseline Interferometric (VLBl) observations of Sgr A* have 
not been able to probe its intrinsic structure, due to the scattering by the interstellar 
electrons (e.g. Davies et al, 1976; Lo et al 1981, 1985, 1993; Backer et al. 1993; Rogers et 
al. 1994; Krichbaum et al. 1997; Bower and Backer 1998). In this paper, we report our 
efforts to image Sgr A* with the Very Long Baseline Array (VLBA) nearly simultaneously 
at five wavelengths (A =6.0, 3.6, 2.0, 1.35 cm and 7 mm). The multi-wavelength imaging, 
with the same interferometer, is crucial for our differentiating interstellar scattering effects 
from the intrinsic source structure of Sgr A*. Finally, more than 20 years after its discovery 
(Balick and Brown 1974), we have for the first time determined that the intrinsic size of 
Sgr A* is asymmetric, being 72 Rsc(*) by < 20 Rsc(*), where Rsc(*) (= 7.5 x 10^^ cm) is 
the Schwarzschild radius of a 2.5x10^ M© black hole. The new VLBA result is consistent 
with the latest millimeter wavelength VLBI results at 3.5 and 1.4 mm, with 3 and 1 
baseline respectively, which suggests that the intrinsic size of Sgr A* may be 17±9 Rsc(*) 
(Krichbaum, T. R etal 1998). 
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2. Observations and Data Analysis 

The observations were carried out using eleven National Radio Astronomy Observatory 
(NRAO)0 25-m radio telescopes: 10 VLBA antennas and 1 VLA antenna. The 7 mm 
observations were carried out on 14 February 1997, while the 1.35 cm observations were 
interlaced with the 2.0 cm observations on 12 February 1997, and the 3.6 cm interlaced with 
the 6.0 cm on 7 February 1997, all at the same hour ranges of UT12h00-20h00. Quasars, 
NRAO 530 and PKS 1921-293, served as amplitude calibrators and fringe detection sources. 
The data recording was in the standard VLBA mode, with 32 MHz bandwidth for both 
circular polarizations at each telescope site. Correlation of the data was done with the 
VLBA correlator in Socorro, New Mexico. All the post-correlation reduction was carried 
out by using the NRAO AlPS and the Caltech VLBl package including Difmap. With 
global fringe fitting, Sgr A* was detected on the short and intermediate baselines depending 
on the observed wavelength, whereas both calibrators (NRAO 530 and PKS 1921-293) 
were detected on all the baselines. The visibility amplitude calibrations were done using 
system temperature measurements at each site. At A7 mm, special attention was paid to 
the elevation dependent opacity corrections, while at A>1.35 cm the atmospheric opacity 
was not significant. Images of Sgr A* were produced for all five wavelengths using standard 
hybrid mapping. 

An elliptical Gaussian model was fitted by the least-squares method to both amplitudes 
and phases in the calibrated visibility data to yield a quantitative description of the source 
structure. Fig. 1 illustrates the model fitting to the 7 mm data, showing both the amplitude 
as a function of visibilities and the closure phase triangles. We emphasize by this figure the 
good quality of the data, the availability of many baseline where the structure can be fitted, 

^The NRAO is operated by Associated Universities Inc., under cooperative agreement 
with the National Science Foundation 
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and the robustness of the fit. The steadily improved performance of the VLB A played a 
pivotal role in the quahty of this data set. The total flux density (S;^), FWHM major axis 
diameter (6'major), FWHM minor axis diameter (6'minor), its axial ratio of minor to major 
axis diameters, and the position angle (P.A.) of the major axis of the model fit are given 
in Table 1, in which we also include previous 7 mm results from 1994.74 for comparison 
(Bower and Backer 1998). At wavelengths >1.35 cm, the mean P.A. of the major axis is 
80° ± 3° (essentially E-W) and the mean axial ratio is 0.53 ± 0.07. At 7 mm, the axial 
ratio, 0.83±0.11, is significantly different from that at the longer wavelengths, while the 
P.A. is not. 



3. The Intrinsic Structure of Sgr A* 

Our near-simultaneous multi-wavelength VLBA mapping of Sgr A* allows us 
to plot both the measured major and minor axis diameters versus the observing 
wavelength in Fig. 2. The measured major axis diameters (open circles) can be fit by 
^major = (1-43 ± 0.02) A^''^^^"^ marc s (A in cm), represented by the sohd fine. Within 
the accuracy of the experiment, such an index is indistinguishable from 2, so that the 
major axis diameters appear to follow a A^ law over a range of A from 7 mm to 6 cm, in 
excellent agreement with the previous result of 1.42 A^ ° (e.g. Alberdi et al. 1993). This 
A^-dependence is consistent with the measured size being dominated by the scattering angle 
that is a result of the radiation from Sgr A* propagating through the interstellar medium 
with fiuctuations in the electron density. The power spectrum of the density fiuctuation of 
interstellar electrons is normally assumed to be oc k~^, where k is the wavenumber of the 
irregularities. The scattering angle scales as X^+'^/(P-'^\ where (1 + 2/{(3 — 2)) = 2, 2.2 for 
/3 = 4, 11/3 respectively (Romani, Narayan, and Blandford 1986). 

Along the minor axis, a fit to the measurements at all five wavelengths yields 



-6- 



^minor = (1.06 ± O.IO) X^-'^^^^-^'^ maic s which appears inconsistent with interstellar 
scattering. A fit to all points for A >1.35 cm, however, yields a dependence of 
(0.87±0.23) X^-^'^^^-^^ marc s, which is consistent, within the errors, with the A^-dependence 
expected for interstellar scattering. If we assume A^-dependence derived for the major axis 
also apphes to the minor axis, ^minor = (0.76 ± 0.05) A^° for A > 1.35 cm (cf Fig. 1). This 
dependence also agrees very well with ^major = (1.43 ± 0.02)A^'^^^°'°^ and a constant axial 
ratio of 0.53. 

This is the first time the X- dependence of the minor axis diameters is determined 
directly by observations, the results of which strongly suggest interstellar scattering 
dominates the observed minor axis image size at A >1.35 cm. The elongation of the 
scatter-broadened image can be caused by an anisotropic scattering medium in the vicinity 
of the Galactic center. The anisotropy in the electron fluctuations in the interstellar medium 
has been postulated to be due to turbulence in a magnetized plasma, and the mechanism 
to generate the density fluctuation has been proposed to be due to the specific entropy 
being mixed by shear Alfvenic turbulence that has "eddies" elongated in the direction of 
the magnetic field on small spatial scales (Higdon 1984; Goldreich and Sridhar 1995). 

Importantly, Fig. 2 shows a significant deviation at 7 mm between the measured 
minor axis diameter, ^minor, and the scattering angle, ^sc, extrapolated from 0.76 A^ °: 

= ^minor - ^sc = (0.58 ± 0.07) marc s -(0.37 ± 0.02) marc s = (0.21 ± 0.07) marc 
s. The deviation of ^^obs from the A^ dependence is naturally expected when the intrinsic 



source diameter, Oi^t: becomes comparable to the scattering angle, since 6'ininor = y^fnt + ^sc 
(Narayan and Hubbard 1988). Thus, this implies that at 7 mm, ^int = (0.45 ± 0.11) marc 
s for Sgr A* along the minor axis (P.A. — —10°; nearly N-S) direction. If we use 
the measurements of Bower and Backer (1998), ^^minor = (0.55 ± 0.11) marc s, 
A9 — (0.18 ± 0.11) marc s and 6'int — (0.41 ± 0.17) marc s. Combining the two sets of 
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measurements, we obtain ^minor = (0.57 ± 0.06) marc s, = (0.20 ± 0.06) marc s and 
^int = (0.44 ± 0.09) marc s. We note that this intrinsic size scale is larger than the value 
inferred at 1.4 mm by Krichbaum etal (1998). In addition to a possible A dependence on 
the intrinsic size, we note that the 1.4 mm measurements were based on a single baseline. 

The reasons that we can discern the intrinsic source size at A 7 mm, contrary to 
previous expectations that intrinsic source size is observable only at A < 1 mm, are (1) 
the scattering size along the minor axis is half that along the major axis, and (2) the near 
simultaneous multi- wavelength mapping of Sgr A*, with the same instrument, over the same 
hour angle, and calibrated in a uniform manner, makes it possible to extrapolate precisely 
the minor axis scattering angle at A 7 mm. 

Conceivably, changes in the refractive properties of the interstellar medium could lead 
to the deviation indicated above, since the 7 mm refractive scattering time scale for Sgr A* 
(proportional to A^) could be short: trcf = O^J^ /Y < 0.5 x (10 kms~VV) year, where D, 
and V are the distance to the scattering medium, the relative velocity of the scattering 
medium and the observer, respectively. However, since the 7 mm source parameters did not 
change over 2.4 years and probably over a longer period, the deviation of the minor axis 
size from 0.76 A^" is unlikely to be due to changing refractive scattering effects at 7 mm. 
Furthermore, other evidence suggests that the refractive scattering effects for the Galactic 
Center must be very small (Romani, Narayan, and Blandford 1986). 

Along the major axis (P. A. of 80°; essentially E-W) direction, the measured diameter 
of 0.7 ±0.01 marc s and the extrapolated scattering size of 0.69 ±0.01 marc s imply that the 
intrinsic size along the same direction has to be < 0.13 marc s. Combined with the minor 
axis intrinsic diameter of 0.44 ± 0.09 marc s derived above, this implies that tie intrinsic 
source structure of Sgr A* could be elongated along an essentially N-S direction, with an 
axial ratio of < 0.3. This also implies an intrinsic brightness temperature of>1.3xlO^°K. 
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Up till now, there exist only limits to the intrinsic size of Sgr A*. Given that the 
intrinsic source size of Sgr A* can now be estimated at 7mm, we can also ask whether there 
are constraints on the wavelength dependence of the intrinsic source size. At A3. 5 mm, the 
upper limit to the observed size is 0.2 marc s, from which we can infer an upper limit to the 
intrinsic size along the minor axis of < 0.18 marc s (Rogers et al. 1994). At A1.4 mm, the 
marginal detection with an interferometer with a fringe spacing of ~0.3 marc s suggests a 
size scale of 0.05-0.15 marc s (Krichbaum et al. 1997; 1998). From the absence of refractive 
scintillation due to focusing and defocusing of the scattered image by large scale plasma 
fluctuations (Gwinn et al. 1991) at A1.3 mm and 0.8 mm, the respective lower hmits to 
the intrinsic size are 0.02 marc s and 0.008 marc s. Taken all together, a A"-dependence of 
^int, with 1.9 > a > 0.7 is not inconsistent with the above limits, at least along the minor 
axis direction. Clearly, this very preliminary determination of the wavelength dependence 
of the intrinsic source size has to be improved with further observations. 

4. The Emission Mechanism for Sgr A* 

There have been several models for the structure and mechanism of radio emission 
from Sgr A*. They typically involve synchrotron emission: from a pulsar wind that is 
confined by ram pressure (Reynolds and McKee 1980), from thermal electrons heated by 
the dissipation of magnetic energy as the mass-loss in the winds from stars in the vicinity 
of Sgr A* such as IRS16 is spherically accreted by a massive black hole (Melia, Jokipii, and 
Narayanan 1992; Melia 1994), from a jet in a coupled jet-disk system (Falcke, Mannheim, 
and Biermann 1993), or from the thermal electrons at an electron temperature of ~ 10^'^ K 
of a two-temperature plasma in a rotating advection-dominated accretion flow (Narayan et 
al. 1998). 

Until now, in the absence of the intrinsic source structure, for confirmation the various 
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models rely on comparison to the spectral energy distribution (SED) of Sgr A* from the 
radio wavelength to the 7-ray (Melia 1994; Narayan et al. 1998; Duschl and Lesch 1994; 
Serabyn et al. 1997; Falcke et al. 1998). However, because of the insufficient angular 
resolution in the wavelength bands shortward of Al mm, it is uncertain that all the radiation 
in the SED actually originates from Sgr A*, making the comparison less than definitive. 
In contrast, the intrinsic source structure derived here provides direct spatial constraints 
on the various models. Specifically, the spherical accretion model (Melia 1994) predicts a 
A 7 mm size too large to be consistent with the results here, while the coupled disk-jet 
model (Falcke, Mannheim and Biermann 1993) predicts a jet size scale smaller than the 
72Rsc(*) obtained here. The advection dominated accretion fiow (ADAF) model for Sgr A* 
can naturally account for the luminosity far below that implied by the estimated accretion 
rate (Narayan et al. 1998). However, if the A 7 mm radiation from Sgr A* originates as 
thermal synchrotron radiation from the inner part of the advection dominated accretion 
disk, the current model has difficulty explaining the elongated shape with an axial ratio 

< 0.3 and the brightness temperature of > 1.3 x 10^° K. An additional component, a radio 
jet, may be needed in the ADAF model to account for the intrinsic elongation of Sgr A*. 
Obviously, observations of the intrinsic source structure at A < 7 mm will probe ever closer 
to the event horizon of the massive black hole and will provide further constraints and 
stimuli for the models of Sgr A*. 

Given the proximity of the center of our Galaxy, understanding the radio emission 
from Sgr A* provides a unique opportunity for probing the physical conditions to within 

< 20Rsc(*) of a 2.5x10^ M© black hole. 

We thank R. Narayan and T. Chiueh for helpful discussions. Z. Shen thanks the Su-Shu 
Huang Astrophysics Research Foundation for travel support to the Array Operations Center 
for data reduction. Research at the ASIAA is funded by the Academia Sinica in Taipei. 
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Fig. 1. — A plot shows the Gaussian eUiptical model fitting (sohd curves) to the visibihty 
data at 7 mm (vertical bars). Left panels: amplitude vs. baseline pairs (BR-FD, BR-KP, 
FD-OV, FD-KP, OV-Y, and FD-LA). The maximum baseline lengths of the pairs are 2346, 
1914, 1508, 744, 1025 and 607 kilometers, respectively. BR-FD is the longest basehne in NS. 
Right panels: closure phase triangles. 

Fig. 2. — A log-log plot of measured (FWHM) source size versus observing wavelength for 
Sgr A* (7- 14 February 1997). The solid line represents a 1.43 X^-^^ fit to the major axis sizes 
(open circles), while the dashed line a 0.76 ° fit to the minor axis sizes (filled circles). 
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Table 1. Pcirameters of Elliptical Gaussian Model Fit 



A 




s. 


^major 


^minor 


Axial 


P.A. 


(cm) 


(GHz) 


(Jy) 


(marc s) 


(marc s) 


Ratio 


(°) 


6.03 


4.97 


0.60±0.09 


49.6±4.50 


25.1±2.00 


0.51±0.09 


81±3 


3.56 


8.41 


0.73±0.10 


18.0±1.53 


9.88±1.68 


0.55±0.14 


78±6 


1.96 


15.3 


0.68±0.06 


5.84±0.48 


3.13±1.14 


0.54±0.21 


73±14 


1.35 


22.2 


0.74±0.04 


2.70±0.15 


1.50±0.59 


0.56±0.25 


81±11 


0.69 


43.2 


1.03±0.01 


0.70±0.01 


0.58±0.07 


0.83±0.11 


87±8 


0.69t 


43.2 


1.28±0.1 


0.76±0.04 


0.55±0.11 


0.73±0.1 


77±7 



^ For comparison, the corresponding results at 43.2 GHz from 1994.74 (Bower and Backer 1998). 
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